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ABSTRACT 

We present evidence for the presence of a transition in the accretion properties of 
radio-loud sources. For a sample of radio galaxies and radio-loud quasars, selected 
based on their extended radio properties, the accretion rate is estimated from the 
black hole mass and nuclear luminosity. The inferred distribution is bimodal, with 
paucity of sources at accretion rates, in Eddington units, of order ~ 10~ 2 - assuming a 
radiative efficiency of 10 per cent - and possibly spanning one-two orders of magnitude. 
Selection biases are unlikely to be responsible for such behavior; we discuss possible 
physical explanations, including a fast transition to low accretion rates, a change 
in the accretion mode/actual accretion rate/radiative efficiency, the lack of stable 
disc solutions at intermediate accretion rates or the inefficiency of the jet formation 
processes in geometrically thin flows. This transition might be analogous to spectral 
states (and jet) transitions in black hole binary systems. 
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1 INTRODUCTION 

Radio-loud active galactic nuclei (AGNs) are characterized 
by the presence of relativistic jets, whose energetics at least 
in some sources is comparable to that produced by the 
nuclear engine. Why and how some systems are able to 
form powerful jets remains one of the unsolved issues in 
AGN research, and currently limits our understanding of 
the relation between the active nucleus, its host galaxy and 
the larger scale environment. It is possible that all AGNs 
pass through a jet phase (evolutionary interpretation) or 
that only certain objects develop powerful jets. For fur- 
ther progress to be made, it is crucial to determine if and 
how the nuclear activity, in terms of accretion properties 
onto the central black hole (BH) and/or the black hole 
spin, is related to the presence and formation of jets. MHD 
processes appear to be most likely involved in the extrac- 
tion (and col limation) of energy to for m jets from the ac- 
cretio n disc l|Bla ndford & , Payn e) Il982l) and/or black hole 
spin felandford fc Znaiekl Il977l) . A significant amount of 
theoretical work has been carried out recently in order to 
understand the effectiveness of these processes - which is 
largely related to the strength and topology of the magnetic 
field - in relation to the properties/regime of accretion (e.g. 
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A further relevant element tightly related to this is- 
sue is the evidence that jets appear in two rather neatly 
distinct flavors, the more powerful ones exhibiting FR II 
radio morphology ; the weaker ones FR I morphology 
jFanaroff fc Rilevlll974<) . At the zeroth order, the unified 
model for radio-loud AGNs identifies these two typologies 
of radio galaxies (RGs) with the parent populations (i.e. 
the misaligned counterparts with respect to the jet direc- 
tion) of radio-loud quasars (RLQs) and BL Lac objects 
respectively, whose observed properties are dominated by 
the highly anisotropic emis sion from relativistic emitting 
jets felandford fc Reeslll97gt) . This association is supported 
by the similar isotropic properties as well as the relative 
number density of such (parent and beamed) populations 
JUrrv fc Padovanilll995l for a review). For the high power 
systems the lack of broad emission lines in FR II RGs has 
been interpreted as due to the presence of material obscuring 
the line of sight at large angles with respect to the jet axis 
(torus). The existence of an obscuring torus has found direct 
support from the detection of scattered broad emission lines 
and from infrared studies which have revealed quasar nuclei 
hidden by rest frame visual extinction in severa l FR lis (e.g. 
lAntonuccilll984l; iHill. Goodrich fc DePovlll996t) . 

Although this scenario is successful at the zeroth or- 
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der, issues are still open. On a rather robust observa- 
tional ground, it has emerged that the fraction of broad 
lined objects increases with the source luminosity, at 
least at low luminosities, and thus - in flux limited sam- 
ples - with redshift (however this appears not to be 
the case when we ak narrow line objects are excluded, 
IWillott et all 12001ft . Furthermore the quasar narrow line 
([OIII]) luminosity and submillim eter emission are lar ger 
than those of radio galaxies (e.g. IWillott et al.l l2002al for 
a recent overview). The two main interpretations of such 
findings invoke either a dependence of the torus open- 
ing angle on source luminosity ('receding t orus model', 
lLawrencdll99ll iHill. Goodrich fc DePovlll996ft or the pres- 
ence of a population of - low luminosity - FR lis with 



low levels of nuclear accretion emission IHinc & Loneaii 
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aspect is the recent suggestion that there might exist a so far 
undetected population of QSOs with FR I radio morphology 
jRlundell fc R.awlingsilzOOll). 



Furthermore, despite of the clue offered by the findings 
bv lLedlow fc Owenl l)l99fift that the radio power at the tran- 
sition between FR Is and FR lis is correlated to the host 
galaxy magnitude, it is still unclear what causes the FR I- 
FR II dichotomy. Possibilities include the different ways in 
which a jet interacts with the am bient medium (iBicknelll 
ll995l:lGopal-Krishna fc Wiital20oJft and/or the different nu- 
clear intrinsic properti es of the accretion and jet forma- 
tion and the jet content l|Baum. Zirbel fc Q'Deall995llMeieil 
Il999l: iRevnolds et al.lll996ft . While little observational evi- 
dence has been gathered in support of any of these inter- 
pretations, the use of correlations between the global large 
scale and nuclear pro perties for the sample analyzed by 
iLedlow fc Owenl l|l99fift intriguingly suggests that a funda- 
mental parameter determining the FR I-FR II dichotomy 
could be the mass accretion rate in Eddington units, which 
possibly regulates the accretio n mode and also the char - 
acteristics of nuclear outflows jGhisellini fc Celottil EoOlft . 
Therefore, while the environment could play a role in de- 
termining the radio morphology and the jet appearance, it 
is of fundamental importance to determine any connection 
between the accretion properties and formation of jets. 

In this context the aim of our work is to shed light on 
the process of formation and evolution of jets with respect 
to the nuclear AGN properties. In particular we intend to 
assess the role of the central black hole mass (Mbh) and the 
mass accretion rate (in Eddington units, rh = M/M^dd) in 
determining the nuclear and extended properties for a sam- 
ple of radio-loud objects (RGs and RLQs) as complete as 
possible. Clearly, m relies on an estimate of M, which will 
be (necessarily) inferred from the bolometric accretion lumi- 
nosity of the AGN, Lboi- The unknown radiative efficiency 
of the accreting flow will constitute a crucial parameter in 
interpreting our findings. 

The outline of the paper is the following. In Section 2 
we introduce the considered samples of sources; in Section 3 
we discuss the analysis and assumptions adopted to ob- 
tain the central black hole mass, the bolometric luminos- 
ity and the mass accretion rate in Eddington units. The 
results and their discussion are presented in Section 4 and 
5, respectively, while in Section 6 we summarize our con- 



clusions. For comparison with previous work in this field 
a zero cosmological constant universe with Ho=75 km s _1 
Mpc -1 and qo=0.5 is adopted [Assuming a ACDM cosmol- 
ogy (J7a = 0.7, £Im = 0.3) does not alter our results, and 
quantities can be easily re-estimated from Table 1]. 



2 SAMPLE SELECTION 

Our goal is to investigate the nuclear properties of a well 
defined sample of radio-loud objects, whose selection is un- 
biased with respect to the nuclear properties themselves. 
Our starting point are the 298 RGs and the 53 RLQs from 
the 3CR catalogue dSpinrad et alJll985ft . 

For each source, we must be able to estimate Mbh (from 
the stellar luminosity of the host galaxy, Mbui ge , see §3.1), 
as well as the bolometric luminosity of the active nucleus. 
While Mbuige is easily measured for RGs, in RLQs the stellar 
light is often overwhelmed by the non-thermal contribution 
of the AGN. From the complete sample of 53 3CR RLQs we 
have therefore selected t he sub-sample for w hich the host 
galaxy is clearly resolved l)Lehnert et alll99fift . This includes 
31 objects with 0.3< z <2.0 and total radio luminosities at 
178 MHz in the range 10 34 < L 1T8 < 10 36 erg s" 1 Hz" 1 . 1 

On the other hand, while Lboi is relatively easy to es- 
timate (modulo an appropriate bolometric correction) for 
RLQs, as the emission from the central AGN is bright and 
easily measurable, the situation is more delicate for RGs. To 
obtain the best possible constraints on Lboi, we select the 79 
3CR radio galaxies for which high resolution Hubble Space 
Telescope (HST) images led to either a measurement or an 
upper limit of the optical unresolved nuclear luminosity 2 
(Chiaberge, Capetti & Celotti 1999, 2000, 2002). This sam- 
ple includes 80% of the 3CR RGs with z < 0.3. Of these, 26 
are FR Is (10 31 < L 17S < 10 35 erg s" 1 Hz" 1 ) and 53 FR lis 
(10 32 < LiTg < 10 35 erg s" 1 Hz" 1 ). While not complete, this 
sample is homogeneous (all galaxies belong to the radio-flux 
limited 3CR sample) and unbiased with respect to nuclear 
properties (see discussion in §5). We will consider this our 
main sample. We also note the significant overlapping of the 
sources both in Lns and z (see Fig. 0. 

In order to assess the incidence of possible selection bi- 
ases, enlarge the statistics, and include RLQs with redshifts 
overlapping with the redshift range of the RGs (thus reduc- 
ing the redshift-luminosity degeneracy and any evolutionary 
effect, see Fig.0, we also consider three additional samples. 
The first comprises the 28 low luminosity RGs (mostly FR I) 
from the B2 s ample (B2RG) with published high resolution 
HST images JCapetti et alJ 12002) . These galaxies populate 
the low redshift range z < 0.15 and span 10 29 < Lns < 10 33 
erg s -1 Hz -1 . The second sample includes the 12 RLQs from 
iMcLure fc Dunlod (1200 ift (MDRLQ). This is a rather het- 
erogeneous sample with redshift 0.15< z <0.3; the common 
thread being the availability of HST images for all objects. 



1 Twelve of the 3CR RLQs were not observed by Lehnert et al. 
because of scheduling problems. 

2 3CR 386 was included in the Chiaberge et al.'s study. However, 
ground based observations subsequently confirmed the putative 
nucleus to be a foreground star (Marchesini, Ferrarese & Celotti, 
in prep). We will exclude the galaxy from further consideration. 
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Figure 1. Total radio luminosity L\-j% at 178 MHz versus red- 
shift. Circles are FR Is (empty ones are B2RGs while filled ones 
are 3CR RGs), squares are FR lis and triangles are RLQs (empty 
ones are 3CR RLQs; light filled ones arc MDRLQs and the dark 
filled ones are the three RLQs from Hutchings et al. 2002). 



The third sample includes t hree high redshift (1.9 < z <2.4) 
RLQs imaged with HST bv lHutchings et all (12002) . 

The total radio luminosity-redshift plane for all of the 
sources considered in this work is shown in Fig. Q 



3 ANALYSIS 

In this section we describe the steps taken to estimate Mbh, 
the AGN bolometric luminosity, and the associated mass 
accretion rate M = Lboi/ec 2 , where e is the radiative effi- 
ciency. 



3.1 Black hole masses 

Only a handful of nearby quiescent galaxies, and only 
four of the RGs and RLQs in our sample, namel y 
3CR 270 (NGC 4261, [ Ferrarese. Ford fc Jaffd Il99d) . 
3CR 272.1 (NGC 4374. iBower et alJ Il998>) . 3CR 274 
(M 87.lMacchetto et al.lll997l) and B2 2116+26 (NGC 7052, 
Ivan der Marel fc van den Boschll99ctl . have Mbh measured 
directly through dynamical modeling of stellar and gas 
kinematics. For these galaxies, Mbh correlates tightly 
with the stellar velocity dispersion of the host bulge and 
(with si gnificantly larger scatter) with the host bulge mag- 
nitud e llKormendy fc Richstonel Il995l : iFerrarese fc Merritd 
l2000l: iGebhardt et al.ll2000l) . Both c orrelations have been 
shown to hold for nea rby AGNs jFerrarese et alJ l200lt 
iMcLure fc Dunlorj f2002). and we will make the fundamen- 
tal assumption that they can be extended to the RGs and 
RLQs considered in this paper. In particular, since measure- 
ments of the stellar velocity dispersion are available for only 
18 of our objects (of which only four are recent measure- 
ments), we will employ the following characterization of the 
Mbh — Mbuige relation: 



log Mbh = -1.58(±2.09) - 0.488(±0.102)M B , 



(1) 



a high quality m easurement of the m ass of the central 
black hole exist s iGeb hardt et al.ll2003t and references in 
iMerritt fc Ferraresdl200ll) . The reduced \ 2 of the relation 
is 4.7, significantly smaller than the one {\ 2 = 23) reported 
for the complete sample of gal axies, including lenticular s 
and spirals, with measured Mbh iFerrarese fc Merrittl2 000). 
Because of worries about the reliability of Mb u i ge measure- 
ments in spirals, and since all of the galaxies considered in 
this paper are genuine ellipticals, we will use eq. 1 for our 
analysis. The scatter in eq. 1 corresponds to a 42 per cent 
uncertainty in Mbh, slightly larger than the value (30 per 
cent) quoted bv lMcLure fc DunlorJ l)2002j) based on a now 
out-of-date sample. 

For the 3CR sources, rest-frame V— band abso- 
lute m agnitudes are taken from the references listed in 
IZirbel fc Bauml l)l995l) or, when not available, from the 
apparent V— band magnitudes reported in the 3CR cata- 
logue. Corrections for Galactic extinction were estimated 
mostly from Bu rstein fc Heiles (1982, 1984) and - w hen not 
possible - from ISchleeel. Finkbeiner fc Davisl <ll998l) . In all 
color correction B — V=0.96, appropriate for e llipti- 
cal g alaxies, was applied llFukueita. Shimasaku fc Ich ikawa 
1995) . For the galaxies in the B2 sample, M B is from 
llmpev fc Greeo rini (1993). i?-band magnitudes of the host 
galaxies of the 3CR RLQs h ave been estimated from HST 
images (ILehnert et al I I1999I) from the residual luminosity 
after subtraction of the nuclear PSF, within an aperture of 
radius l".4 (which, at the redshift of the targets, covers the 
entire galaxy). For the RLQs in the McLure & Dunlop and 
Hutchings et al. samples, Mr are listed by the authors. A 
color correction B — R=1.57 (Fukugita et al. 1995) was ap- 
plied to transform R— to B— band magnitudes. Cosmologi- 
cal and K corrections have been applied to all magnitudes. 

For each galaxy, the rest-frame absolute magnitudes 
and the BH masses (from eq. 1) are listed in Table 1. 



where Mbh is in solar masses and M B is the total B— band 
luminosity of the host elliptical galaxy. This is derived for 
the most current sample of 16 elliptical galaxies for which 



3.2 Bolometric luminosities 

All of the objects in our sample have HST images. With 
the exception of the RLQ sample of McLure & Dunlop, the 
optical nuclear luminosity, L , is derived by fitting the unre- 
solved nuclear component with the appropri ate HST point 
spread function (Chiaberge et al. 1999, 2002; ICapetti et alJ 
l2002l:lHutchings et alJ2002t) . For the RLQ sample of McLure 
& Dunlop, L is obtained fr om the rest-frame monoc hro- 
matic luminosity at 5100 A jMcLure fc DunlodlioOlh . K- 
corrections have been applied in all cases. Uncertainties on 
L vary for the different samples, but are typically of the 
order of 15 per cent. 

Estimating the nuclear bolometric luminosity Lboi from 
L involves the knowled ge of the bol ometric correction, rj, 
which we take from Elv is et alJ (119941) . This was derived for 
a sample of radio-loud and radio-quiet quasars, and might 
not be appropriate for the radio galaxies considered here. We 
will discuss the possible consequences of our assumptions in 
§5.1. 

Nuclear bolometric luminosities are listed in Table 1 
for all objects. The uncertainty on the adopted r\ is ~35 per 
cent, yielding an average error for Lboi of ~ 38 per cent. 
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4 RESULTS 

In the following, we will discuss whether different classes of 
objects in our sample have different properties in terms of 
BH mass, nuclear bolometric luminosity and/or mass accre- 
tion rate. For radio galaxies, we will make the distinction 
between FR I and FR II. Furthermore, in view of the evi- 
dence that the FR II population might be heterogeneous in 
terms of its nuclear properties (see §5) the FR II RGs will be 
further divided according to their optical spectral properties 
in low- and high-ionization emission line galaxies (LEGs and 
HEGs, respectively) and broad line RGs (BLRGs). 



4.1 Black hole masses 

Although the objects in our sample span about three decades 
in Mbh (Fig. most BH masses are in the 10 s - 1O 9 M 
range. The paucity of objects (at least radio galaxies) host- 
ing BH with smaller mass is possibly due to selection ef- 
fects (althoug h the samples were sele cted bas ed on ra- 
dio power, see iMcLure fc Dunlodl2002t but also lHoll2002t 

m 



lOshlack. Webster fc Whitindl2002l : IWoo fc Urrvll2002alh 

The upper panel in Fig.|2]shows the distribution of black 
hole masses Mbh for all 3CR sources (FR I, FR II and RLQ); 
the lower panel shows the distribution of Mbh for all of the 
objects considered in this study. 

A Kolmogorov-Smirnov statistical test (Table 2) shows 
significant differences in the distribution of Mbh only be- 
tween RLQs and both FR Is and FR lis, but with mean 
values within their standard deviations. No evidence for dif- 
ferences are found for objects belonging to different sub- 
classes or different samples (3CR RGs vs B2 RGs, or 3CR 
RLQs vs RLQs from McLure & Dunlop). 

4.2 Bolometric luminosities 

Fig. [3] shows the distribution of nuclear bolometric lumi- 
nosity Lboi for FR Is, FR lis and RLQs in the 3CR sample 
(upper panel) and for all of the sources in this study (middle 
panel). Unlike the distribution in Mbh, objects belonging to 
different classes show systematic differences in their nuclear 
(optical or bolometric) luminosity (see Table 2). This is nei- 
ther a new nor an unexpected result (indeed the distinction 
between RGs and RLQs is based on their nuclear luminosity, 
see also §5). However, it is intriguing to notice that the over- 
all distribution - spanning about seven orders of magnitude - 
appears structured, with two peaks easily recognizable, and 
a shortage of sources at Lboi ~ 10 44 — 10 45 erg s _1 . The 
low luminosity peak (Lboi ~ 10 42 ' 7 erg s _1 ) is comprised in 
equal measure of FR I and FR II sources, while the high 
luminosity peak includes mainly radio-loud quasars. While 
there are no FR I galaxies straddling the two peaks, some 
FR lis extend to the high luminosity regime populated by 
the quasars. 

We estimated the statistical signi ficance of a possible bi- 
modality using the KMM algorithm i Ashman. Bird fc ZepJ 
Il994f) . The distributions for the entire sample, and for the 
3CR sample alone (both RGs and RLQs) are both strongly 
inconsistent with being unimodal (P-value <0.05, see Ta- 
ble 2). The distribution of the 3CR radio galaxies alone, 
however, is consistent with unimodality (P-value >0.1), al- 
though there is marginal evidence for a second peak around 




Figure 2. Distributions of Mbh for the 3CR sources (upper 
panel) and for all of the sources considered in this work (lower 
panel): the thin continuous line is the distribution for all of the 
objects, the horizontal shaded area represents the FR Is, the ver- 
tical one represents the distribution for the RLQs and the thick 
continuous line indicates the FR lis. 



10 44 erg s _1 . This peak is populated exclusively by FR II 
sources; if these are further divided into LEGs, HEGs and 
BLRGs, an additional trend emerges: all LEGs and about 
2/3 of the HEGs are low luminosity sources, while 1/3 of 
the HEGs and all of the BLRGs have high nuclear luminos- 
ity (see Fig. |31 lower panel). We will discuss this point in 
more detail in §5. We point out that sources for which only 
upper limits to the nuclear luminosity exist, are found ex- 
clusively in the low luminosity peak (see Table 2), thus not 
affecting these findings (and the equivalent ones about the 
rh distributions). 
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Figure 3. Distributions of Lboi for the FR Is, FR lis and RLQs 
sources belonging to the 3CR sample (upper panel) and for all of 
the objects considered in this work (middle panel): the thin con- 
tinuous line is the distribution for all of the objects; the horizontal 
shaded area represents the FR Is, the vertical one represents the 
RLQs and the thick line refers to the FR lis. The lower panel 
shows the distributions of Lboi f° r the LEG (horizontal shaded 
area), HEG (thick continuous line) and BLRG (vertical shaded 
area) sub-samples; the thin continuous line is the distribution for 
all of the objects. 



4.3 Mass accretion rates 

The mass accretion rate in Eddington units is then easily 
estimated as m = Lboi/e^Edd- The resulting distribution is 
shown in Fig.0]for the 3CR sample (upper panel) and for all 
the considered objects (middle panel), under the assumption 
(to which we will return in §5) that the radiative efficiency 
e is equal to unity for all sources. 

A KMM test (Table 2) shows that the overall distri- 
bution in rh is non unimodal (P- value <0.05), a clear con- 
sequence of the fact that, while the distribution in Mbh is 
unimodal, the one in Lboi is not. In this case, the rh distri- 
bution for the 3CR radio galaxies alone is also inconsistent 
with being unimodal, (P-value <0.05; see Table 2), although 
the Lboi distribution for the same sample is consistent with 
being unimodal. Note that the typical uncertainty on rh 
(Alogm ~0.45) is smaller than the width of the peaks and 
significantly smaller than their separation. All of the BLRGs 
and 1/3 of the HEGs have rh >2xl0~ 3 , while all LEGs and 
2/3 of the HEGs have rh <2xl0" 3 . The LEGs have sta- 
tistically the same distribution in rh as the FR I galaxies 
(see Table 2), while HEGs have an asymmetric distribution, 
peaking at rh ~3x 10~ 4 and rapidly declining toward smaller 
rh. 

As a final test on the presence of bimodality we also 
considered the effect of the dispersion in the Mbh — Mb cor- 
relation from which we estimated the BH masses. By adopt- 
ing a Gaussian probability distribution of Mbh (for a given 
Mb) with width equal to the scatter of such relation we 
found that the bimodality is still significant, but only when 
all of the 3CR sample (and not only the radio galaxies) is 
considered. 



5 DISCUSSION 

The main findings of this work are: a) within the sample 
considered, FR Is, FR lis and radio-loud quasars are indis- 
tinguishable based on A/bh which power the nuclear activity, 
however b) they not only differ significantly in both nuclear 
bolometric luminosity and mass accretion rate (as several 
previous works have already pointed out), but there is evi- 
dence for a bimodal behavior in such properties. We stress 
that the errors associated with Lboi and Mbh translate, as 
mentioned before, to an uncertainty on the mass accretion 
rate smaller than the width of the peaks. It is conceivable 
that samples with smaller uncertainties on Mbh and Lboi 
would produce an even more pronounced multi-modal dis- 
tribution in rh than the one shown in Fig. [I] 

As shown by Chiaberge et al. (2002) , the nuclear optical 
core luminosity appears to be markedly different for FR I 
and most of FR II radio galaxies. The low FR I optical lu- 
minosity well correlates with the radio one, suggesting that 
the emission in the two bands has a similar origin, namely 
non-thermal radiation from the jet. Furthermore, the high 
fraction of detected nuclear sources imply the absence of a 
geometrically thick obscuring torus, as also supported by 
the lower level of mi d and far IR emission and limited UV 
absorp tion in FRIs (iHeckman et alJll994t | Chia berge et alJ 
2002b|). FR II sources show instead a more complex behav- 
ior, as expected given the already proposed view - based on 
the optical spectral properties - that the FR II class is not 
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Figure 4. Distributions of rh for the FR Is, FR lis and RLQs 
belonging to the 3CR sample (upper panel), for all of the sources 
studied in this work (middle panel): the thin continuous line rep- 
resents all the sources, the horizontal shaded area indicates the 
FR Is, the vertical shaded area the RLQs and the thick line repre- 
sents the distribution for the FR lis. The lower panel is the same 
as the middle one, but only LEGs and BLRGs are represented by 
the thick continuous line, while HEGs are excluded. The reported 
values correspond to a radiative efficiency e = 1. 



a homogeneous population llLaing et al ESI). Some FR lis 
follow the same correlation as FR Is, while others show a 
significant 'excess' in their optical luminosity, given their 
radio power. This excess has been interpreted as the contri- 
bution from a thermal component asso ciated with radiation 
from accretion (see also the results bv lWillott et alj |2002b 
on the anti-correlation between the 4000 A Ca break and 
the line luminosity). This hypothesis is strongly supported 
by the fact that the different dominance of non-thermal vs 
thermal contribution correlates with the optical line prop- 
erties: LEGs behave like non-thermal FR Is, while BLRGs 
show a thermal excess. HEGs show both behaviors. As sug- 
gested by the large line equivalent width of HEGs with ap- 
parent FR I— like optical luminosity, and in agreement with 
the basic unification scenario, this is consistent with HEGs 
being thermal BLRG-like objects whose optical emission is 
partly obscured (Chiaberge et al. 2002). In this scenario, 
low-ionization FR II and FR I sources belong to the same 
population, while high-ionization FR lis and BLRGs cor- 
respo nd to misaligned radio-loud quasars llJackson fc Walll 
1999). This unification scheme also accounts for the find- 
ings that some BL Lacs have both extended radio mor- 
phology dKolkaard et al]ll992l) and extended radio power 
ilCassaroer^nTl999r) typical of FR II sources, as well as 
for the behavior of the luminosity function and evolution of 
radio galaxies (e.g. iDunlop fc Peacock! Il99d : IWillott et alJ 
Il999t Grimes et al. 2003). 

How do our findings fit into this overall picture? In the 
above scenario, the dichotomy between non-thermal (jet) 
and thermal (accretion) dominated sources is reflected into 
the observed bimodal distribution in Lboi- This implies that 
all of the HEGs would shift toward higher Lboi if their in- 
trinsic nuclear luminosity could be observed: although the 
amount of absorption is not known a priori, the above unifi- 
cation scenario implies that the HEG intrinsic luminosity is 
comparable to that of BLRGs/RLQs (as appears to be the 
case for the couple of sources with spectropolarimetric data, 
Chiaberge et al. 2002, see also §5.1). Furthermore, because 
the observed optical luminosity of FR Is and LEGs is domi- 
nated by the non-thermal component, it formally represents 
an upper limit to the thermal accretion luminosity of these 
sources. If the two contributions could be separated, and 
only the thermal accretion optical luminosity were plotted 
in Fig. [3] all low luminosity sources would shift toward even 
lower I/boi 3 (see however below), thus exacerbating the 'sep- 
aration' between the two populations of thermal and non- 
thermal dominated objects. 

Before proceeding any further, it is essential to analyze 
possible observational biases which might be responsible for 
the observed behavior. We begin by pointing out that only 
two of the 31 3CR RLQs are flat spectrum. Although for 
these objects Lboi might be contaminated by non-thermal 
emission, removing them would not affect our findings. Fur- 
thermore, observational biases are expected to play against 
the brightest and more distant RLQs, rather than the weaker 
ones, since resolving the underlying galaxy was a necessary 
condition for selecting these objects. It is therefore unlikely 
that our sample is missing low luminosity RLQs which would 



3 Note also that some of the nuclear luminosities in the low- 
luminosity range already represent upper limits, see Table 1. 
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populate the transition regime in either Lboi or m. On the 
other hand, no bias is expected against bright optical radio 
galaxies in the complete 3CR sample observed with HST, 
since the selection of these sources is based on the extended 
radio power, which does not correlate well with the nuclear 
optical luminosity. Furthermore, the distribution of the ther- 
mal and non-thermal sources in terms of radio power is con- 
tinuous: the two classes overlap, and show no sign of a 'gap'. 

The observed bimodality in the distribution of m should 
also not be induced by possible cosmic evolution effects for 
the following reasons. First, the distribution in m is bimodal 
also when the sample of 3CR RGs alone is considered; since 
these sources are all within z <0.3, and since the FR I and 
the FR If sub-samples overlap both in redshift and extended 
radio luminosity distributions, evolutionary effects (if as- 
sumed to be similar to those affecting RLQs) are not crucial 
for the 3CR RG sample. Second, the lower and the higher 
redshift samples of RLQs have similar properties (in). Fi- 
nally, even the cosmic evolution of the high-z 3CR RLQs, 
estimated by using the parameterization L op t oc (1 + z) 3 
l|Bovle fc Terlevicr]IT99sFl between z = 1 and z = 0.1, is a 
factor ~6 (and ~20 for the few z ~ 2 sources) , much smaller 
than the separation between the two peaks of the distribu- 
tion in m and therefore not enough for being responsible for 
the observed bimodality. 

A further possible s ource of bias could be present if - as 
suggested bv lLawrencd (Il99lf) (see also Simpson 1998, 2003) 
- there is an inverse dependence between the degree of ob- 
scuration and the source luminosity ('receding torus model') 
within the FR1I/RLQ population. This would imply a natu- 
ral decrease in the number of detectable sources at lower rh 
and on average a lower luminosity for the obscured objects. 
While there is some observational support to this scenario 
(see the above papers), there are some indications that the 
bias thus introduced would not alter our results. We first no- 
tice that our sample is selected on the basis of the extended 
radio properties, and while indeed some of the sources have 
not been looked at by HST, this was not related to their 
optical nuclear properties - in fact the unobserved objects 
include both HEGs, LEGs and BLRGs in similar number 
(only a few for each sub-class). Secondly, as discussed be- 
low, the bolometric luminosity estimated for HEGs from 
their line emission is compatible with that of BLRGs, as also 
supported by inferences for objects with spectropolarimetric 
data and by the distribution of the core dominance param- 
eter of BLRGs, HEGs and LEGs (Chiaberge et al. 2002). 
This indicates that sources for which line information are 
available (i.e. more than half of HEGs) are not intrinsically 
less luminous. Although these arguments certainly cannot 
rule out a luminosity dependence of obscuration, we notice 
that in any case the expected average luminosity correction 
between radio galaxies and RLQs - to account for the se- 
lection effects induced by the 'receding torus' - is typically 
estimated of order ~ few (compatible with the obs erved ra- 
tio in line luminosity of the two populations, see ISimpsonl 
120031 : Fradhunter et alJll993l) . We conclude that such an ef- 
fect could not be the primary cause of the found bimodality 
nor can significantly affect our results. 

To summarize, selection effects should not qualitatively 
change the observed distributions, and in particular the lack 
of sources at intermediate Lboi- 



5.1 Bolometric corrections and Lboi 

We consider the assumptions made in deriving Lboi next. 
The value of the bolometric correction, 77, adopted relies 
on an average spectral energy distribution for RLQs. On 
the basis of the unification model, it is plausible to assume 
that this bolometric correction is appropriate for BLRGs 
and HEGs as well, which are believed to be intrinsically 
identical to RLQs. However, while there are indirect indica- 
tions for low accretion luminosities in low power RGs and 
BL Lac objects (from the lack of strong emission lines in 
these systems, from limits on the amount of nuclear photon 
densi ties and from the detected level of X-ray emissi on, e.g. 
ICelotti. Fabian fc Reedll99St iDi Matteo et al]|2000h . there 
is poor observational and theoretical information on their 
spectral energy distribution and even on the nature of the 
emission of the accretion component. Therefore it is not at 
all obvious that a bolometric correction which is valid for 
RLQs can also be applied to FR I and low excitation FR II 
radio galaxies. 

However, there are two additional ways of estimating 
Lboi- 

The bolometric 'disc' luminosity can be estimated 
from the observed line luminosity, assumed to originate 
from gas photoionized by the accretion radiation field. 
This approach has several drawbacks, including the fact 
that both the shape of the underlying photoionizing 
continuum, and the covering factor of the line emitting 
gas are unknown in these systems. A comparison between 
Lboi estimated via a bolometric correction, and through 
the observed line luminosity can be carried out for 49 
radio galaxies with accurate and homogeneous (in quality 
and available lines) [Oil] and/or [OIII] line luminosity 
information (see the on-line database maintained by 



Chris Willott, 
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. Assuming that the line radia- 



tion follows gas photoionization by the nuclear continuum, 
from the line luminosities we estimated the total luminosity 
in narrow lines (as Lnlr = [3(3L[ n] + 1.5L[ nr])], e.g. 
IRawlines fc Saunders] Il99 lT) . and in turn the ionizing 
luminosity Lion = LnlrC -1 , where C is the covering 
factor. For a typical C ~ 0.01 for all the LEGs (with 
the exclusion of 3C 123, which could be a mis-classified 
HEG, Chiaberge et al. 2002) and the BLRGs we find a 
good agreement - within a factor of a few - between the 
estimates of Lj on and Lboi (inferred from the assumed rf). 
Consistently with the discussion above, HEGs have instead 
Li on comparable with that of BLRGs (and thus largely 
in excess of their Lboi). Seven (out of eleven) FR Is show 
also good agreement between the two estimates, while the 
remaining four reveal an excess of Li on by a factor 30-300. 
In order to further assess this issue, we also chec ked (via 
modeling using the code CLOUDY. lFerlandlll99a) on the 
possible difference in the narrow line ratios (oxygen lines 
vs total luminosity) which could arise if the (unknown) 
spectrum of FR Is was significantly different from the 
typical quasar-like spectrum. We find that even in the cases 
of spectra as hard as F(y) oc u° the above estimates of 
Lion are robust (within factors of order unity). Significant 
deviations are only achieved for harder spectral shapes or 
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values of the ionization parameter so low (U <-4.5) to be 
insufficient to give raise to the observed [0111] luminosities. 
While the reason for the observed discrepancy for the 
four FR Is is not clear at this time, we notice that for 
all four objects Lnlr has been derived from the [Oil] 
line lumino sity only (by assuming a r atio of [OII]/[OIII] 
fluxes, see iRawlings fc Saundersl Il99ll) which could be 
more significantly contaminated by emission not associated 
to nuclear p hotoionization (such as shocks induced by 
the jets, e.g. |Popitalll995l or the line contribution from 
the galaxy itself, see iBaum. Zirbel fc O'Deal I1995T) . Note 
that the similarity of the luminosities inferred by the two 
approaches in low power objects also suggest that (at 
least for those sources) the levels of non-thermal and 'disc' 
emission are in fact of the same order, reinforcing the 
presence of low luminosity/low rh peaks rather than tails in 
these distributions. 

For ten sources (9 RLQs and 1 BLRGs) we also esti- 
mated Li on f rom the luminosity in broad lines (fol lowing the 
procedure in|Cdotti^Pad£vani&, Ghisellinilll997t line lumi- 
nosities from |j^ick^on^^B^wn^^9lTr iinding again agree- 
ment with Lboi within a factor 1.5. Given the uncertainties 
inherent to both methods, we conclude that there is a good 
agreement between the luminosities Lboi estimated via the 
bolometric correction rj, and those estimated via the line 
emission. 

We should, however, point out that estimates of nu- 
clear luminositie s based on line emission fo r a sample of 
BL Lac objects tWang. Stauber fc HollioOSfl - which have 
typical extended radio luminosities similar to our sample of 
RGs - are typically a factor 30 larger than those of our RGs 
based on the assumption of a 'universal' bolometric correc- 
tion. This discrepancy is pos sibly due to the fact that while 
Wan g. Stauber fc Hoi i2002T) used both narrow and broad 
emission lines to estimate the bolometric luminosity, they 
adopted a covering factor typical for broad line emission 
(C = 0.1). 

It is also relevant to notice that no clear gap has 
been fo und in the line emiss i on properties of radio g alax- 
ies (e.g.lZirbel fc Baumlll995t IBaum. Zirbel fc 0'Dealll995l 
IWillott et aLl Tl999) where the low frequency radio emission 
has been found to correlate with the narrow line luminosity. 
Nevertheless the results by the above authors clearly show 
the presence of a transition (change of slope in such corre- 
lation) between FR Is and FR lis at values comparable to 
those inferred here, which indicate that the line luminosity 
(for the same radio one) in FR lis is larger by a factor > 
10 res pect to what observed in FR Is. In fact IWillott et all 
(2000) discuss the possibility that an analysis including the 
spectroscopic information on the line excitation might lead 
to a clearer separa tion. Furthe r sup port to this comes from 
the recent work bv I Wills et alJ (120041) . who targeted a sam- 
ple of 13 low-luminosity RGs. The correlation between nar- 
row line [OIII] luminosity and radio power in FR II galaxies 
has the same slope, but significantly different zero point, 
than the one observed for FR I galaxies. In particular, the 
latter have systematically lower optical line luminosity at a 
gi ven radio power ; in ag reement with the results obtained 
bv lZirbel fc Bauml dl995l) and Baum et al. (1995). 

The second, unfortunately even more uncertain, way 
of estimating Lboi is through theoretical arguments. If the 
accretion regime in FR Is and LEGs corresponds to an op- 



tically thick, geometrically thin, radiatively efficient flow, 
the expected spectral energy distribution could be approxi- 
mated at first order as a multicolor blackbody, plausibly sim- 
ilar in shape to that of RLQs, but with a temperature and 
flux scaled to reflect the much lower mass accretion rates. 
In this case, the appropriate bolometric correction would 
be a factor ~ 4 lower than the value derived for RLQs, in 
the sense of moving the lower end of the distribution to 
lower bolometric luminosities, thus making the separation 
between non-thermal and thermal sources even more pro- 
nounced. However, as amply discussed in the literature, at 
low accretion rates the accreting flow could attain a two tem- 
perature, optically thin, geometrically thick flow structure. 
Although the dominant physical processes, and the resulting 
conditions of the flow in such a regime are still unclear, emis- 
sion due to cyclo-synchrotron, Compton and bremsstrahlung 
emis sion by a thermal electron pop ulation could be expected 
(e.g. lNaravan. Barret fc McClint ock|l99 7|). For instance, the 
most developed spectral model llQuataert fc Naravanlll999l) 
predicts a bolometric correction which is about a factor ~ 
6 larger than the one used in our analysis (for a typical set 
of parameters and m ~ 7 x 10 ) due to the contribution 
of the low energy cyclo-synchrotron component. This does 
work in the sense of 'filling the gap' between non-thermal 
and thermal sources, i.e. partly washing out the bimodality 
in Lboi- [Note that the (partial) agreement of the luminosi- 
ties estimated via line emission and Lboi does not affect this 
possibility, as the major component of the luminosity in this 
case would not contribute to the gas photoionization]. We 
however stress that this possibility inherently assumes that 
the apparent bimodality arises as a consequence of a change 
in the flow radiative efficiency. Indeed, as we will discuss in 
§5.2 the reduced efficiency associated with an optically thin 
flow obviously acts in the sense of 'closing the gap' in the rh 
as well as in the Lboi distributions. 

We conclude that in terms of their bolometric luminos- 
ity, radio-loud sources are naturally divided in BLRG/RLQ, 
which show a relatively broad peaked distribution at the 
high Lboi end, and FR Is plus low-excitation FR lis, pop- 
ulating the low Lboi peak. Selection effects do not seem to 
account for the deficit of sources at intermediate Lboi- The 
observed distribution can be reconciled with an underlying 
unimodal distribution apparently only if the bolometric cor- 
rection for the low Lboi sources has been underestimated by 
a factor of several, as could be the case if e.g. accretion in 
these sources proceeds through a geometrically thick, opti- 
cally thin accretion flow. 

The fact that the BH mass does not play a role in deter- 
mining the different nuclear optical properties of the radio- 
loud populations, the observation that the distribution in rh 
show signs of being more strongly bimodal than the one in 
Lboi, as well as the fact that rh is ultimately regulating the 
emitted luminosity, suggest that rh can be the key parame- 
ter responsible for the systematic difference in the emission 
of thermal and non-thermal sources. 

5.2 Radiative efficiency vs rh 

We now discuss the role of the mass accretion rate in driv- 
ing the bimodal behavior of the different populations. It is 
important to stress that the sources with Lboi < 3 x 10 44 
erg s _1 - i.e. in between the two peaks of the rh distribu- 
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tion - are HEGs, i.e. objects partly absorbed, which should 
be intrinsically as luminous as BLRGs. This would imply 
the presence of an underlying gap in the distribution span- 
ning about one/two decades in luminosity (see Fig. [3] lower 
panel) and in accretion rate (few 10 -4 < rh < fewlO -2 - see 
the lower panel of Fig. 0] in which HEGs are excluded). The 
presence of such a wide separation between the two peaks 
would thus require a physical process able to account not 
only for a transition, but also for a significant gap in the 
observed quantities. 

A first possibility is that rh is actually discontinuous, 
in the sense that for m below a certain value (rh ~ 0.01) 
the evolution time-scale of accretion is much faster than 
for higher rh. This would account for a fast decrease in 
the number density of 'thermally' dominated sources and 
a peak/tail at much lower rh. This could correspond to 
an on-off switching of the nucl ear accretion between dif- 
ferent regimes in rh (e.g. also iDi Matteo et al] l2003t) or 
an evolution ary sequence from the FR II to the FR I- 
like regime JZirbel fc Bauml Il995t iLedlow fc Owenl Il996l 
iGhisellini fc Celottill200lt ICavaliere fc D'Eliall2002D . 

Alternative possibilities where all m are possible, 
equally probable and have a continuous distribution are 
however open. The transition could then be due to a 
change (but not simply a smooth decline as predicted by 
simple advection flow models) in the radiative efficiency e 
of the accreting flow. Theoretically, the transition between 
a radiative efficient, geometrically thin, optically thick 
accretion flow, and a radiatively inefficient, geometrically 
thick and op tically thin configuration is expected t o occur 
at m ~ a 2 \Rees et alJll982t INaravan fc Yilll995lh For a 
viscosity parameter a ~ 0.1, considered to be phenomeno- 
logically and physically plausible (e.g. iMahadevanl Il997t 
iBalbus fc HawlevHl998th this value is indeed close to the 
transition in the rh distribution shown in Fig. |1] in this 
scenario the peak at low rh is an artificial consequence 
of our assumption that e is the same for all sources, 
while in fact the radiative efficiency for the objects at low 
m is much smaller (and possibly a decreasing function 
of rh itself, INaravan fc Yil Il995t) . However our limited 
sample and the theoretical uncertainties do not allow to 
expand further on this point, nor (yet) to give any robust 
indication on the valu e of a. A simila r interpretation 
has been suggested by 1 Wills et al] (|2004) in relation to 
the offset in the [OIII] luminosity between FR Is and 
FR lis. However, another plausible possibility in a two 
temperature optically thin configuration is that the flow 
is subject to significant energy /m ass loss, in the form of 
an ou tflow, as p ointed out bv | Bla ndford fc Begelmarj 
J1999D (see also lHawlev fc Balbusl 120021: and for a 
different view iNaraxam Garcia fc McClintockl Il997t 
INaravan . leumeiishdie ^&^^bramowicS\^00^tl or a jet 
(e.g. iPellegrini et all I2003D . In such" a case in fact the 
actual rh accreting onto the BH would be much smaller 
than that at large disc radii. Note also that a bimodal 
distribution in m, and especially the presence of a gap, 
might conflict with the proposed possibility of an increase 
- if continuous - of magnetic dissipation of disc energy into 
the corona (possibly associated wi th the formation of jet) 
at decreasing mass accretion rates llMerloni fc Fabianll2002t 
iLivio. Pringle fc KindEoolj) . Clearly in these situations, as 



discussed in §5.1, also the intrinsic distribution in Lboi can 
of course be continuous. 

Different scenarios can also be envisioned, as at values 
of m ~ 10 -2 other physical processes could be relevant. For 
instance, in the optically thick, high surface density disc so- 
lutions a transition from the gas pressure to radiation pres- 
sure dominated accretion regime for supermassive BH can 
occur, for 'reasonable' values of the viscosity parameter: this 
could thus imply the presence of an un stable flow in the in- 
ner d isc onsetting at these values of m ( Sh akura fc Sunvaevl 
Il973l) . While a paucity of sources over a narrow range in rh 
would disfavor such possibility, an actual gap could in fact 
be consistent with this view. 

Alternatively, the intermediate region could correspond 
to a range of rh for which the formation of a powerful jet 
is inhibited. In MHD models of relativistic j et for mation, 
and in particular in the Bla ndford fc Znaiekl il977t) one in 
which the black hole rotational energy is extracted via a 
Penrose-like process, the key parameter determining the ef- 
ficiency of the process is the stren gth of the poloidal mag- 
netic field Bp . If, as stressed by ILivio. Qgilvie fc Pringlel 
(1999), a realistic condition for the equilibrium poloidal 
magnetic field is B p ~ (H/R)Bj, (where H/R is the in- 
ner height to distance ratio in the disc and B$ is the az- 
imuthal component of the magnetic field), then the process 
of jet formation could be efficient enough only for H/R ~1, 
i.e. either for rh <Cl (accretion-starved advective disc) or 
for rh >1 (su per- Eddington advective-disc), as suggested by 
M eicr ( 2002). Again, while the presence of a significant num- 
ber of sources at intermediate rh may cast doubts on such 
interpretation, the assumption that HEGs are intrinsically 
identical to BLRGs, leads to an absence of sources in a range 
of almost two decades in rh. 

As a final note, we point out that for e ~ 0.1 some RLQs 
result to be super-Eddington, although by only a factor of 
a few. In order to ascertain whether these RLQs could be 
partly affected by relativistic beaming, we estimated their 
radio spectral indices, and found that only five (out of the 17 
RLQs above the Eddington value) are actually flat spectrum 
sources, indicating that relativistic beaming likely does not 
largely account for the highest luminosities. Also, the ma- 
jority of these sources are not the highest redshift RLQs in 
our sample (six are at z <0.3 and only five are at z >0.8), 
and thus evolutionary effects do not largely account for their 
extreme luminosities. 



6 CONCLUSIONS 

We considered the nuclear properties of radio-loud objects, 
comprising a complete sample of radio galaxies (both FR I 
and FR II) and radio-loud quasars. We estimated the mass 
of the central supermassive black hole, the nuclear accretion 
luminosities, and the corresponding mass accretion rates. 

We found that while BH masses span a relatively large 
range, they are clustered around 10 8 — 10 9 Mq and there 
is no systematic difference in their values between FR I, 
FR II radio galaxies and r adio-loud quasars, in agreemen t 
with previous findings (e.g. lO'Dowd. Urrv fc Scarpall2002T) . 
Biases introduced by the (radio) criteria adopted in defin- 
ing the samples do not allow us to address whether a 
lower limit on the BH mass or an upper envelope in ra- 
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dio luminosity for a Riven mass exist in radio-loud sou rces 
jMcLure fe DunlorJ200llLaorl2000UDunlop et alJ2003t but 
see iHol 12001 lOshlack. Webster fc Whiting! 12003: Woo & 
Urry 2002a, 2002b). 

However, the distributions in bolometric luminosity, 
and (even more so) in the corresponding mass accretion rate 
appear to be bimodal. FR I, low excitation FR II and some 
high excitation FR II galaxies occupy the low end of the 
Lboi and rh distributions: 10 41 < Lboi < 10 43 ' 5 erg s _1 
or 10 -6 < rh < few 10 -4 . On the other hand, radio-loud 
quasars, broad line radio galaxies and some high excita- 
tion FR II galaxies occupy the high end of the distribu- 
tions: 10 45 < Lboi < 10 47 erg s" 1 or 10 -2 < rh < 1. The 
region in between is characterized by a marked deficiency 
of sources, especially if high excitation galaxies are indeed 
intrinsically similar to broad line ones. The relative com- 
pleteness and homogeneity of our sample, which includes 
low and high power sources selected on the basis of the 
extended radio properties, allowed to reveal such bimodal 
behavior, which had not been noticed in previous studies 
jO'Dowd. Urrv & Scarnal EEl IWoo fc TIrrvl l2002ah. This 
is so pronounced to appear despite of the uncertainties in 
the inferred quantities. No clear selection bias which could 
induce such result has been found. In a future work, we 
will test this bimodal behavior through more accurate de- 
termination of the BH masses (via the Mbh — cr relation) 
and luminosities for a smaller subsample of sources includ- 
ing a representative sample of FR II sources (Marchesini, 
Ferrarese & Celotti, in preparation). The range of rh corre- 
sponding to the transition/gap between FR Is and FR lis 
(rh ~ 10 -2 — 10 _1 for e = 0.1) is similar to the value of 
rh ~0.06 corresponding to the dividin g line between FR I 
and F R II radio galaxies discussed bv lOhisellini fc Celottil 
teOOll) . which was inferred via an independent proce dure 
based on the properties of the lLedlow fc Owenl il996f) RG 
sample via global correlations between the extended radio 
emission and the accretion one. In particular the magni tudes 
of the host galaxies reported bv lLedlow fc Owenl jl996T) were 
converted into BH masses (again via the Afb u i gc — Mbh rela- 
tion), while the radio luminosities were associated to nuclear 
accretion rates, by adopting the c orrelation between radio 
power and narrow line luminosity (Wil lott et al.l fl999L and 
the latter converted into i onization luminosity. The a bsence 
of a gap in the findings bv lGhisellini_j£_Cel_ottj| j 200j|Vmight 
however arise as the relation bv IWillott^tHd" ll 19991) was 
(necessarily) adopted by these authors also for FR Is, which 
tends in fact to overestimate t he FR I line luminosity (as 
also shown bv I Wills et ai1l2004t) . 

A dichotomy is not seen if the radio luminosity is sub- 
stituted to Lboi (in this case the distribution appears contin- 
uous with significant overlap between the different classes), 
making our findings even more intriguing. A discontinuity in 
the rh distribution implies a transition/gap in the accretion 
which is not reflected in the extended radio emission and 
possibly, if a similar jet radiative efficiency would character- 
ize FR Is and FR lis, in the jet power. This would also im- 
ply a significantly different ratio between the accretion and 
jet powers in low and high luminosity systems, questioning 
whether the rh taken at face value in low power sources is 
sufficient to give raise to the associated jets. We also notice 
that the radio luminosity corresponding to the transition 
in L op t between LEGs and BLRGs - according to the es- 



tablished radio-line correlations (e.g. IWillott et aljfl999h - 
refle ct the break in the radi o sources luminosity function 
fe.g. lDunlop fc Peacocklll990D . 

Our findings indicate that the whole radio-loud popula- 
tion can be characterized by thermal vs non-thermal dom- 
inated (optical) emission (as suggested by Chiaberge et al. 

2002) . The population of FR II radio galaxies is in this re- 
spect complex and - as expected - the source behavior is 
related to the nuclear line properties. In other words, the 
nuclear structure and the central activity define the parent 
radio galaxies of RLQs and BL Lacs and such activity is not 
uniquely connected to the extended radio power and mor- 
phology. This might suggest different time-s cales for the evo- 
lution of the radio and optical activity f see ICirasuolo et alJ 

2003) and/or a temporal evolution from the higher to the 
lower power sources. 

Some possibilities accounting for the found bimodal be- 
havior have been discussed. These involve a relatively rapid 
transition in the accretion rate; a transition in the accretion 
mode, implying a change in radiative efficiency and/or the 
onset of significant mass loss, with transition values of rh in 
agreement with the rates at which radiative inefficient flows 
are expected to be possible; the presence of an unstable (e.g. 
radiation pressure dominated) accretion regime or an ineffi- 
cient jet formation process associated to geometrically thin 
flows at intermediate rh. While no definite conclusion can 
be drawn on the origin of the found transition/gap, some of 
these hypothesis might be (dis)-favored by examining with 
more accurate measures the effective absence of sources over 
a significant range at intermediate rh, i.e. the presence of a 
significant gap in the observed properties. A further aspect 
which we are exploring is whether the dichotomic behavior 
holds also for radio-quiet sources, which would imply that 
the last of the above possibilities cannot be the (only) ex- 
planation. 

Several pieces of evidence already indicated that in 
low power radio galaxies and BL Lac objects the lumi- 
nosity resulting from dissipation in the accretion flow is 
extremely low, compared to that typically inferred from 
RLQ activity (as e.g. from the l ow line luminosity, lim- 
its on the nuclear emission, e.g. lUrrv fc Padovanil [ l995|: 



Celotti. Fabian fc Reesil998l:IO'Dowd. Urrv fc Scarpal2002t 

Wang. Stauber fc HolbOO^ Furthermore, stringent limits 



came with the finding that the radiated output is or- 
ders of magnitude lower than that expected from the 
(Bondi) accretion rate estimated from the properties of 
hot gas at the accretion rad ius for a 10 per cent ef- 
ficien cy (iFabian fc Reesl I199E ; e.g. M87, iDi Matteo et alJ 
l2003t IC4296. iPellegrini et al.ll200^ . Our findings, consis- 
tent with the above results, appear thus to indicate the 
presence of a relatively rapid transition and a gap between 
the high radiative efficient accretion regime and the ac- 
cretion properties of low power sources. Clearly, these re- 
sults are also co n sisten t with the possibility, suggested by 
IDi Matteo et alJ a2003T> . that the low disc accretion rate 
might have a time dependent (on-off) character (even be- 
tween low activity and almost 'inactive' states). Note that if 
indeed a large fraction of the radio, X-ray and optical emis- 
sion in some low power objects is dominated by non-thermal 
radiation from the base of the jet (as proposed), this might 
imply rh lower than ~ 10 -5 — 10~ 6 . The findings from the 
line emission however indicate that at least in a relevant 
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number of sources the non-thermal jet and 'disc' emission 
should be of similar order. 

Finally, it is worth mentioning that the found behavior 
could indeed be analogous to the transitions between spec- 
tral states present in black hole X-ray binary systems and 
microquasars, possibly in terms of the value of the transi- 
tional m, the nuclear spectral pr operties and t he onset and 
nature of the jet production (e.g. lFenderll2003T) . 

As already amply stressed in the literature in relation 
to the understanding of accretion at low rates, there is great 
need to test for the presence of reflection features in the 
X-ray spectra of such systems, as well as variability and 
polarization meas ures, as it has already been possible for th e 
Galactic Center feaganoff et al.l l200ll : lAitken et alJl200(J) . 
Theoretical and numerical work might however lead to a 
deeper understanding first. 
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Table 1. Properties 


of the sources 
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46.74 



Table 2. Results of statistical tests 



Kolmogorov-Smirnov Test p 



Pop. 1 


Pop. 2 


M B H 


Lbol 


m 


FR I 


FR II 




7xl0 -8 


2xl0 -7 


FR I 


RLQ 


2xl0" 3 


3xl0~ 22 


3xl0~ 22 


FR II 


RLQ 


5xl0- 7 


5xl0- 20 


lxlO" 17 


3CR FR I 


B2RG 








LEG 


HEG 








LEG 


BLRG 




7xl0~ 6 


7xlO~ B 


HEG 


BLRG 




5xl0~ 8 


7xl0~ 7 


3CR RLQ 


MDRLQ 








FR I 


LEG 








FR I 


HEG 




lxl0~ 4 


3xl0~ 5 


FR I 


BLRG 




5X10" 10 


5xl0~ 10 


3CR FR I 


3CR FR II 




3xl0~ 4 


lxl0~ 4 






Means and standard deviations 


Population 


log M B h 


log Lbol 


log rh 




FR I 


8.7+0.4 


42.30+0.83 


-4.51+0.79 


B2RG 


8.8+0.3 


42.16+0.81 


-4.72+0.75 




FR II 


8.6+0.5 


43.53+1.09 


-3.14+1.03 




LEG 


8.5+0.7 


42.70+0.64 


-3.89+0.69 




HEG 


8.5+0.4 


43.19+0.82 


-3.41+0.91 


BLRG 


8.7+0.5 


44.92+0.50 


-1.91+0.30 




RLQ 


9.1+0.5 


46.13+0.60 


-1.11+0.36 








KKM Test 




Parameter 


Ml 




a 2 


P-value 






All radio-loud sources 




log Lbol 


42.60 


45.77 


0.76 


<0.001 


log m 


-4.30 


-1.49 


0.54 


<0.001 






3CR radio-loud sources 




log Lbol 


42.85 


45.76 


0.79 


<0.001 


log rh 


-4.09 


-1.60 


0.49 


<0.001 






3CR RGs 




log Lbol 


42.58 


44.42 


0.51 


0.101 


log rh 


-4.20 


-2.39 


0.43 


0.032 



Mbh and Lbol are in units of M n and erg s 1 , respectively. Ac- 
cording to the KMM algorithm iAshman. Bird fc Zep3ll994tl m 
and (j.2 are the means of the two Gaussian distributions represent- 
ing the bimodality, while <r 2 is their common covariance; when the 
P-valuc is <0.05, the observed distribution is strongly inconsis- 
tent with being unimodal; when P-value G [0.05,0.1] the observed 
distribution is only marginally inconsistent with being unimodal; 
therefore, when P-value is >0.1, the observed distribution is not 
inconsistent with unimodality. For the Kolmogorov-Smirnov sta- 
tistical test, only values for the significantly different populations 
are reported. 
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